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We analyze the full counting statistics (FCS) of quantum dots in the Kondo regime contacted 
by normal and superconducting leads or an STM tip. To describe the Kondo resonance we use 
an effective model for the quantum dot in the Kondo regime in combination with the FCS for 
tunneling contacts calculated using the Hamiltonian approach. We show that the situation of 
weak coupling to the superconducting electrode in the case of superconductor gap being larger 
than the Kondo resonance width can be easily handled and verify the method by comparing our 
theoretical predictions to experimental data. This allows us to make predictions for the noise and 
cross-correlation in setups involving a superconductor. We find that a positive cross-correlation is 
possible in the case of a superconductor and two normal leads contacted via two quantum dots in 
the Kondo regime. 
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I. INTRODUCTION 

Accessing electronic transport in nanoscale systems re- 
quires an understanding not only of the current but also 
its fluctuations. Whereas charge transfer between normal 
metals proceeds via single electron transfer processes, 
in normal-superconductor hybrid structures the current 
transport at energies below the superconducting gap is 
governed by Andreev reflection^ An electron from the 
normal conductor cannot enter the superconductor at en- 
ergies below the gap and is therefore retroreflected as a 
hole leaving a Cooper pair behind. 

In recent years it has become possible to not only 
make nanometer sized junctions 2,3 but also to contact 
a single impurity via two leads.-' 5 This configuration 
and the corresponding Anderson impurity model allows 
to study the Kondo effect— in many geometries^— 
Normal-quantum dot-superconductor hybrid structures 
are of special interest since they represent the main build- 
ing blocks of Andreev entajiglersj 1 

To obtain a complete knowledge of the charge transfer 
in these systems it is not sufficient to study the nonlin- 
ear I — V and the shot noise alone. It is very useful to 
possess a complete information about the probability dis- 
tribution of charge transfer, also referred to as full count- 
ing statistics (FCS). It is usually calculated in terms of a 
counting field A-dependent cumulant generating function 
(CGF) of charge transfer i 19 : 20 Derivatives of the CGF 
with respect to the counting field directly generate all ir- 
reducible moments of the distribution. Additionally the 
type of distribution and the actual A-dependence allows 
for a rigorous identification of the elementary processes 
of charge transfer Although the experimental access 
to higher order cumulants is rather difficult it has already 
been demonstrated in a number of works i 25 ' 26 

The FCS for quantum point contacts (QPCs) be- 
tween normal and/or superconducting leads have been 



considered before . ' — Likewise the FCS for a quan- 
tum dot in the Kondo limit contacted by normal drains 
have been calculated at the Toulouse pointi 32 ' 33 How- 
ever, these techniques are hard to generalise to the case 
of a superconducting drain. Likewise the analytical tech- 
niques developed so far for the I — V characteristics of 
normal-quantum dot-superconductor device a 11 : 34 even in 
the simplest cases lead to rather cumbersome calcula- 
tions. 

In this work we want to show how results for the 
FCS of tunneling junctions allow for effective descrip- 
tions for the situation of a quantum dot in the Kondo 
limit that is coupled to a superconductor or an STM 
tip. We analyse the interplay of the Kondo effect 5 - and 
s up er conduct ivit y 1 1 ' 36 for situations in which Tk < An, 
where An is the superconductor gap and Tk is the Kondo 
resonance width. We also assume a strong asymme- 
try of the couplings of the quantum dot to the normal 
electrode and the superconductor/STM tip. While an 
STM tip represents a special normal contact the asym- 
metry is generic in the case of a normal-quantum dot- 
superconductor junction in this regime/ 3 ^ As a first step 
we derive the FCS for a normal-superconductor junction 
using the Hamiltonian approach as introduced in [37j ■ I n 
this approach to quantum transport nanostructures are 
analysed starting from a microscopic Hamiltonian, which 
essentially describes transport as tunnelling contrary to 
the alternative strategies of [28ll29j modelling transport 
by scattering. Our approach allows us to formulate the 
effective model for a coupled quantum dot. We find 
that the explicit consideration of Andreev reflection in 
normal-quantum dot-superconductor junctions improves 
the fit to the experimental conductance compared to pre- 
vious studies and leads to a considerable enhancement of 
the Fano factor for voltages below the gap. We also show 
how Andreev reflections in a junction with two quantum 
dots lead to a positive cross correlation of the currents in 
the normal leads. 



2 



II. QUANTUM POINT CONTACTS 

The first system we want to consider is a QPC between 
a normal metal and a superconductor /STM tip as in 
[29l | . It is described by the following Hamiltonian 



H = H 1 +H 2 +H T 



(1) 



Hi refers to the normal metal, which is modelled as a 
noninteracting fermionic continuum at chemical potential 
fix, written in terms of electron field operators ^kia- H 2 
either refers to another normal metal or a superconduct- 
ing reservoir that is described by the BCS Hamiltonian 38 
with the gap Aq of the superconducting terminal 



H 2 = X>**iL*fc2<r 
A;, it 

k2i + ^-k2l^k2t 



As in previous treatments of the problem^ the second 
lead is kept at /i 2 = 0, so that the applied bias voltage is 
given by V — M2 = Mi ( we use units e = H = ks = 1)- 
Electron hopping between the two electrodes is described 
by the Hamiltonian 39 

H T = ^7 [*+ (1 = 0)tf 2ff (z = 0) + h.c] , 

a 

where 7 is the amplitude of the respective tunneling cou- 
pling and the tunneling is assumed to occur at x = in 
the coordinate system of the electrodes. 

The density of states (DOS) in the vicinity of the Fermi 
level in the normal electrode poi is assumed to be only 
weakly energy-dependent. For a superconducting drain 
the DOS is strongly energy-dependent and diverges for 
energies close to the superconducting gap. From BCS 
theory 3 ^ it is known to be given by 



P2 = P02 



A 2 
^0 



In this work we want to use the Hamiltonian approach, 
meaning the complete resummation of perturbation the- 
ory in the tunneling coupling. Instead of coupling the 
counting field directly to the current operator as in 
[28||29| . A enters the tunneling part of the Hamiltonian 
according to the prescriptio n 24 ! 33 

TX = £ 7 [e~ iX/2 *t(x = 0)* 2ff (x = 0) + h.c] . (2) 



The counting field is explicitly time-dependent: it is 
zero outside the measurement interval < t < r and 
carries different signs on the forward/backward Keldysh 
branches. The CGF is given by the expectation value: 
ln%(A) = ln(e lAQ ) that allows to compute the cur- 
rent, noise and also higher statistical moments that 
may become experimentally observable in these systems. 
The different cumulants of the distribution are obtained 
by successive differentiation with respect to the count- 
ing field. The CGF for long measurement times r 
can be calculated analytically using a Green's function 
metho d 19 ' 40 specifically designed for studying quantum 
impurity problems! 33 ' 41 ' 42 

For the FCS calculation we follow the approach out- 
lined in [33| : for long measurement times the CGF may 
be written using an adiabatic potential^ 3 - (that by con- 
struction does not depend on time) as 



lnx(A) = -i / dtU(X) = -irU(X) . 
Jo 



In turn the adiabatic potential is related to the counting 
field derivative of T x in Eq. © 




£<T C *+ (a: - 0)* 2a (x = 0)) A + h.c. (3) 



where (-)a is defined as (-)\ := l/x(A)(-) with (•) being 
the ordinary expectation value with respect to the sys- 
tem's Hamiltonian in Eq. ([TJ with Ht replaced by T x 
in Eq. The mixed and A-dependent Green's func- 

tions in Eq. ([3]) are calculated exactly by summing up all 
orders in the tunneling coupling generated by T x . The 
FCS for the case of a contact between two normal metals 
is well known^ [see also Eq. (jHJ)]. In the case of a con- 
tact between a normal metal and a superconductor con- 
trary to [33| , of course, we have to use superconducting 
Green's functions for the superconducting side as given 
e.g. in [iij . In this way we obtain an exact expression 
for the temperature and energy-dependent CGF that is 
valid over the whole range of possible parameters. It is 
given by 



lnx7V5(A,r) = t 

+T A2 ( W )(2n 2 - 1) {(2n 2 - 1) [{e lX - 1)V_(1 - n 1+ ) - 2(e lA - l)( e -* A - l)m_n 1+ 
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+ (e- tA - 1)V+(1 - ni-)] + 2n 2 (e* x - l)( e - tA - l)(n 1+ - 1 + m_)} 
+T BC (u;)(2n 2 - l)(e lA - e- lA ) 2 {(2n 2 - l)K_e 4A + n 1+e - tX + T e n 2 (l - n 2 )(e iX - e^) 2 

-(m+ - 1 + m_)n 2 (e a + e- lA )] - 4n 2 (l - n 2 )(n 1+ - 1 + m_)} (^^) } 

+ In {1 + T A {u) [m+(l - m_)(e aiA - 1) + m_(l - n 1+ )( e - 2iA - 1)] } f^r^) } , (4) 



where the effective transmission coefficients are given by 



T e (w) 
TA2(w) 



4r R 



[(i 



r e ) a -r^]' 
[(l + r^-rip 



Tbc{u) 



and 



r* +21^(1 -ra) + (i + raja' 



(5) 



The energy-dependent DOS of the supercond ucting ter- 
minal affects the hybridisations T e = T\uj\/ \J\uj 2 — Aq| 

and T A = rAo/VIA 2 , -w 2 |, where T = 7r 2 p iPo27 2 /2- 
The subscripts e and A refer to the electronic and 
Cooper pair density of states. T e refers to the trans- 
fer of single electrons and Tai, Tbc describe the addi- 
tional noise contributions by Andreev reflection above 
the gap and branch crossing processes respectively— Ta 
is the transmission coefficient of Andreev reflection pro- 
cesses below the gap. The Fermi distribution for the 
normal/superconducting lead is abbreviated by rti + , 2 - 
m_ = 1 — n\ + {— uS) refers to hole-like contributions. 

The expression in Eq. ([J| is the first result of this 
work. For Ao — > only the first line of the expression 
remains which is the sum of electrons and holes of 
the the expression for the CGF of a normal tunnel 
contact derived in 19]. Consequently the dominating 
charge transfer events above the gap are single electron 
transfers. Below the gap only terms proportional to e 2lA 
occur referring to double electron transfers in Andreev 
reflection processes^ Since the result in Eq. Q is 
an exact expression this interpretation is unambiguous 
contrary to perturbative results where the analysis is 
often not triviaL 33 ' 45 Using the scattering approach the 
FCS of a similar system has been studied before i 28 i 29 
Because a direct comparison of the CGF is quite difficult 
we start by comparing the lowest cumulants and find 
excellent agreement. 

In particular, we compute the current corresponding 
to this CGF via 



Ins = — — lnxAr S (A, r)| A=0 - 
t d\ 



(6) 



Considering the differential conductance dl^s/dV for 
low voltages and T = we obtain the result previously 
obtained in [46j . Likewise the equivalence of the widely 



used BTK model 4 ^ and the Hamiltonian approach used 
here has been demonstrated in [37| ■ We also reproduce 
exactly the results for the nonlinear current voltage 
characteristics in [37]. To further demonstrate the 
validity of the Hamiltonian approach we also performed 
a comparison to the experimental results for Al/Cu 
contacts with high transparency presented in and 
obtained a very good agreement. 

We also compared our result directly to the FCS ob- 
tained in [2811 29} . In [28| it was argued that one should use 
the transmission coefficients from the Blonder-Tinkham- 
Klapwijk (BTK) model 4 ^ and in Ref. [2!| circuit theory 
in combination with the BCS Green's functions repro- 
duces this result. In [371 ] it was shown how to obtain the 
BTK transmission coefficients from the result using the 
Hamiltonian approach. Using this mapping on Eq. (f?]) 
we recover the results obtained in [28|29]. From Eq. (f4]) 
it is evident that the FCS has two regimes: for |U| 3> Ao 
only single-electron processes constitute the current while 
for |V| <C A only Andreev processes contribute. 



III. QUANTUM DOTS IN THE KONDO 
REGIME CONTACTED BY 
SUPERCONDUCTORS OR STM TIPS 

The Hamiltonian approach can be applied in more 
complex geometries, e. g. for the resonant level system, 
which is the simplest quantum dot. Including Coulomb 
interaction on the quantum dot we observe a Kondo reso- 
nance when the number of electrons on the quantum dot 
is odd so that it obtains a localized magnetic moment of 
spin 1/2. At energy scales below the Kondo temperature 
Tk the dot spin hybridizes with the lead spin density and 
forms a singlet state. At this strong-coupling fixed point 
a perfectly transmitting channel opens up and the Kondo 
effect can be described as a pure resonant level system as 
far as the electronic transport is concerned^ For weak 
coupling of the second lead to the quantum dot one can 
see from the calculations in [33ll49| that a resonant level 
leads to a renormalization of the effective transmission 
coefficient for a tunnel contact with the density of states 
of a resonant level. This situation emerges in the case 
of a quantum dot in the Kondo regime contacted by a 
normal lead and an STM tip (KNQN) or by a supercon- 
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ductor provided T K < A holds 36 (KNQS). In the op- 
posite case Tk > A the Kondo resonance also strongly 
couples to the quasiparticles on the superconducting side 
leaving the current voltage characteristics almost identi- 
cal to those obtained for normal conducting systems^ 6 - If 
only one drain strongly couples to the quantum dot the 
coupling to the second lead can be described by an ef- 
fective tunneling transmission. We can thus employ the 
transmission coefficients for the QPCs multiplied with 
the Kondo density of states T 2 K / [(u - V) 2 + T 2 K ] M 
If the Kondo impurity is contacted by two normal drains 
with strongly asymmetric coupling its transmission can 
be effectively written as 



4L 



NN 



V 2 
1 K 



(i + r NN ) 2 [(u-v) 2 + r 2 K y 



(7) 



where T P . 



Tattv) 2 , where T 



NN 



2 " 2 is the transmission coefficient of the ordinary 



k/[(«> 



Vf 



r 2 1 

1 K\ 



accounts for the 



4T^/(1 

7^"Pol^oo27 , ' 
QPC. The factor V 2 
Kondo resonance 
Additionally we have to include the density of states out- 
side the Kondo peak given e.g. by the Hubbard subbands 
at ±U/2, where U refers to the strength of Coulomb re- 
pulsion on the quantum dot. In a first approximation this 
can be modelled by a constant background conductance 
Tfc as in (36|. The FCS of these additional processes can 
be described by the standard Levitov-Lesovik formula^ 

ln Xb (A,r) =2r J ^ln{l + T b 

x [(e lA - l)m+(l - n 2 ) + (e- lX - l)n 2 (l - n 1+ )]\&) 

The total CGF is given by \tlxkn = Itixknqn + ba-Xb, 
where xknqn can be derived from \b by replacing TJ, by 

T en K- 

As expected, for the differential conductance we observe 
the typical Lorentzian shape that is given by the form of 
the Kondo resonance. However, in this case there is no 
energy scale but the width of the Kondo resonance. 
This is different in the case of a superconducting lead 
with its characteristic energy gap. E.g. for two super- 
conducting terminals the Kondo effect is suppressed if the 
superconducting gap is larger than Tk j*™ A similar effect 
occurs in the regime IV i$ A considered here in a KNQS 
junction since only electrons on the normal side strongly 
couple to the emerging Kondo resonance. One therefore 
obtains a strong asymmetry of the Kondo coupling be- 
tween the normal and the superconducting drain to the 
doti^ 6 - Consequently, one observes a strong suppression 
of the total conductance while side peaks at the super- 
conducting gap appear. The assembly of superconductor- 
quantum dot-hybrid structures with good coupling to the 
quantum dot makes T k < Aq a typical experimental 



case 
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Since the coupling of the superconductor to the 
Kondo impurity is small and branch crossing and 
Andreev reflection above the gap are both processes 



of higher order we find that we may safely ne- 
glect their contributions. We can thus employ the 
transmission coefficients for the normal-superconductor 
QPC multiplied with the Kondo density of states 
T 2 K / \(u> — V) 2 + T K ~\ for single electron transmission 
and Y K / { [(w - V) 2 + T 2 K ] [(w + V) 2 + I*.] } for An- 
dreev reflections to obtain 



T e K{u) 



4r e 



r 2 

1 K 



(i + r e ) 2 {u-v) 2 + t 



K 



v\ + 2r 2 4 (i-r 2 ) + (i + r 2 ) 2 
x [(u-v) 2 + r 2 K } K [(u + v) 2 + r 2 K ] 



The total CGF, including the background DOS, is given 
by \nxKS = hi-XKNQS + hix&, where xknqs can be ob- 
tained from In xns m Eq. |4] by replacing T e by T e K and 
T A by T A k while setting T A 2 = T BC = 0. 
We calculate the ratio of the differential conductances 
in the case of a KNQS junction and a KNQN device. 
Fig. Q] shows the comparison of the theory outlined above 
with experimental data, where the Kondo conductance 
peak in the KNQN state has been approximated by a 
Lorentzian^ 6 - We observe a fairly good agreement up to 
a small difference for small voltages. This can be at- 
tributed to additional diffusive transport channels. 




FIG. 1: Theoretical prediction for the differential conduc- 
tance Gs-qd-n divided by the normal state value Gn-qd-n 
shown as the solid curve. The parameters taken from [3y] in- 
clude Y K = 0.7A, T = 0.05, T D = 0.053A, T b = 0.375 and 
the width of the Lorentzian for the conductance in the KNQN 
case of 0.6A. To obtain the peaks at the right position for 
our fit we need to employ a slightly lower superconducting 
gap compared to the fit parameter A in the semiconductor 
model^ (dotted curve) taken also from [3(| A = 0.7A. The 
result can be compared to the experimental data shown as 
the dashed curve.— 

In contrast to the prediction by the semiconductor 
model^ 2 - we obtain quantitative agreement for voltages 
up to approximately 2A due to the explicit inclusion of 
Andreev reflection and using the full energy dependence 
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FIG. 2: Theoretical prediction for the Fano factor in the 
normal-quantum dot-superconductor junction with the same 
experimental parameters as used in the fit in Fig. [1] The 
solid curve has been obtained with Andreev reflections and 
the dotted curve is the prediction without Andreev reflec- 
tion (Tak = 0). One observes a clear difference between the 
two curves that amounts up to a difference of about 15% at 
V/A = 0.6 (see dotted line). 

of the transmission coefficients derived in the Hamilto- 
nian approach. The energy-dependent transmission co- 
efficients lead to broader conductance peaks close to the 
superconductor gap while the explicit inclusion of An- 
dreev reflection leads to an enhanced conductance at low 
bias. For higher voltages the structure in the background 
density of states starts to play a role and deviations are 



expected. By introducing a transmission coefficient Tak 
we also included the possibility of Andreev reflection. 
Due to the low transparency of the contact it is strongly 
suppressed and causes only a small contribution to the 
differential conductance at low voltages. However, ac- 
cording to our effective theory they should be observable 
by a noise measurement . 

In Fig. [5] the Fano factor including Andreev reflections 
has been compared to the prediction without Andreev 
reflections (Tak — 0). One observes a clear difference in 
the Fano factors. This effect can be used to identify the 
presence of Andreev processes even for large interaction 
strengths. 

The possibility of Andreev reflection also raises the ques- 
tion whether one could be able to observe a pos itive cross 
correlation in a three terminal setup like in [16M49ll53M54j . 
From the results obtained here we can address the case of 
a superconductor coupled to two Kondo dots contacted 
by two normal leads at the same chemical potential (see 
Fig. p. 

In principle a positive cross correlation of the currents in 
the two normal leads could be inferred from the process 
of crossed Andreev reflection where a Cooper pair from 
the superconductor splits up and the two electrons are 
transferred to the separate normal terminals. The cur- 
rent originating from this process has been analyzed in 
a similar geometry in [55| . For simplicity here we only 
want to address the case of T = 0, where the CGF adopts 
the form 



ln XsNN (X,r) =2rf -ln{l+ £ T eKf} (u)(e^ - 1)}6 ( M_zM +T f ^ ln{1+ £ T AK ^ - 1) 

JO T ^1,2 ^ / J-V 77 /3=12 

+TcAK(e^ + ^-l)}e^^y (9) 
where the transmission coefficients can be inferred using the result from [49] to be 

4r e 



T 



eK/3 



T 



AK0 



TcAK 



(i + r e )* (w-v)3 + iV 

4ri 

v\ + 21* (i - r>) + (i + r 2 ) 2 [( U - vy + r K1 ][(" + v? + r 2 ^] 



8T A r^r 2 



A' 2 



v\ + 21^(1 - yd + (i + yiy [( W - vy + I* X ][( W + Vf + T K2 ] 

I 



The cross correlation can be calculated as a second 
derivative of the CGF 



p i ^ _1 d 2 \u X snn(X,t) 

12 T d\!d\ 2 



(10) 



A=0 



Fig. 2] shows the result for three different configura- 
tions of the widths of the respective Kondo resonances. 



In all cases we observe a positive cross-correlation that 
first increases towards the superconducting gap and then 
quickly dissapears for voltages above the gap. This 
means that in principle it should be possible to observe 
a positive cross-correlation of the currents in the nor- 
mal leads that is mediated by crossed Andreev reflection. 
Above the superconducting gap single electron transfer 
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FIG. 3: Sketch of the experimental setup for measuring the 
cross-correlation. A superconductor (SC) is contacted via two 
quantum dots in the Kondo regime to two normal leads (N) 
at the same chemical potential. 




FIG. 4: Theoretical predicition for the cross correlation P/ 2 
at T — in a three terminal setup with one superconducting 
and a normal lead. All curves have been obtained for V — 0.05 
and a broadening of the BCS density of states has been taken 
into account by a Dynes parameter of F_d = 0.005Ao. The 
solid curve refers to Fki = 0.7Ao = Tk2, the dashed curve is 
for Ta'i = 0.3Ao, Yk2 = Ao and the dotted curve has been 
obtained for Fki = 0.3Ao = Yk2- 



to the separate terminals sets in and causes a negative 
cross-correlation. However, we should mind that the pre- 
diction has been done for T = 0. At finite temperature 
thermally activated single-electron transfer becomes pos- 
sible. This process causes a negative correlation of the 
two currents that will, in most cases, destroy the positive 
cross-correlation by crossed Andreev reflection. Never- 
theless, this shows that even in the Kondo regime a pos- 
itive cross-correlation could in principle be observable. 



IV. CONCLUSION 

We have calculated the FCS for a quantum dot in the 
Kondo regime coupled to a normal drain and a super- 
conductor/ STM tip. We have demonstrated the validity 
of our approach by comparing our results with existing 
studies as well as experimental data. For the normal- 
quantum dot-superconductor junction we found that ac- 
cording to our effective theory Andreev reflections could 
be seen in the noise, leading to a considerable change in 
the respective Fano factor. Finally, we have investigated 
crossed Andreev reflection in a three-terminal setup. We 
have shown that a positive cross-correlation could in prin- 
ciple be observed even in the Kondo situation. 
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